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The Interacting Gluon Model (IGM) is a tool designed 
to study energy flow, especially stopping and leading particle 
spectra, in high energy hadronic collisions. In this model, 
valence quarks fly through and the gluon clouds of the hadrons 
interact strongly both in the soft and in the semihard regime. 
Developing this picture we arrive at a simple description of 
energy loss, given in terms of few parameters, which accounts 
for a wide variety of experimental data. This text is a survey 
of our main results and predictions. 



I. INTRODUCTION 

A. Why a model? 

After more than 30 years of continuous advances we 
might expect that Quantum Chromodynamics (QCD), 
the established theory of strong interactions, would pro- 
vide us with a satisfactory understanding of high energy 
hadronic reactions. Unfortunately this is not yet the 
case and we have to study these reactions using mod- 
els instead of the theory. This is so essentially because 
of two reasons. The first one is because we can only 
perform reliable calculations in the perturbative regime, 
i.e., in reactions where the momentum tranfer is larger 
than a few GeV. However these represent only a small 
fraction of the hadronic cross sections. Indeed, even at 
very large energies most of events involve low momentum 
transfer, as indicated by the average transverse momen- 
tum of the produced particles, which is, in most of the 
experiments, of the order of 1 GeV or less. The second 
reason is that the number of interacting particles may 
be so large that many body techniques and approxima- 
tions are needed. At RHIC, for example, the number of 
finally produced hadrons may be as large as 6000 (!) re- 
sulting from the complicated interaction among a similar 
number of quarks and gluons at the initial stage of the 
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collisions. The study of these systems can not be made 
from first principles and models are required. 



B. What is a good model? 

Since making models is inevitable and since there has 
been a proliferation of models for particle production, we 
must try to establish criteria to decide when a model is 
better than other. A condition to be satisfied by a model 
is a clear connection to the underlying theory, i.e., the 
use of the appropriate degrees of freedom with the cor- 
rect QCD interactions. Moreover, assumptions should 
be made only where the theory is not applicable and the 
introduction of parameters should be restricted to a min- 
imum. Furthermore, a model must have predictive power 
and be testable. Even with these constraints there are 
many implementations of the basic QCD concepts and 
many different ways to treat the non-perturbative dy- 
namics. 

Among all the existing models, there are some which 
try to give a very comprehensive description of all possi- 
ble experimental data. Usually these models are at some 
stage transformed into event generators and are used by 
experimental groups. While they may be helpful in pro- 
jecting detectors and analyzing data, they have the dis- 
advantage of containing many parameters and of being a 
kind of "black box". Well known examples of this type 
of model are HIJING [1], VENUS [2] and NEXUS [3]. A 
comprehensive list of available models of this kind can 
be found in [4] . In a different approach there are models 
which concentrate only on some more specific features 
both of theory and experiment. These models are more 
transparent, easy to handle, have only a small number 
of parameters and are devised to test only some lim- 
ited aspects of the theory. A famous example of this 
kind of model is the thermal model, which is one of the 
first successful models formulated at a very early stage 
of the research on multiparticle production [5] and which 
remains very popular still in our days [6]. This model 
does not involve amplitudes nor cross sections and is ap- 
plied only to systems which might be in thermal equilib- 
rium. Typically we fit simple thermal distributions to the 
experimentally measured transverse momentum spectra 
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and extract the effective temperature. Here we have lit- 
tle input and little output but we may learn something 
studying different systems and, for example, establish the 
behavior of the temperature as a function of the colli- 
sion energy. The model discussed here belongs to this 
cathcgory of "economic" models. Our aim is to decribe 
energy flow (stopping, energy deposition and leading par- 
ticle spectra) with a simple picture based on QCD, with 
few parameters and learn something from the analysis of 
data. 



C. Why study energy flow? 

Multiparticle production processes are the most com- 
plicated phenomena as far as the number of finally in- 
volved degrees of freedom is concerned. They also com- 
prise the bulk of all inelastic collisions and therefore are 
very important - if not per se then as a possible back- 
ground to some other, more specialized reactions mea- 
sured at high energy collisions. The large number of de- 
grees of freedom calls inevitably for some kind of statis- 
tical or hydrodynamical descrition when addressing such 
processes. All corresponding models have to be supple- 
mented with information about the fraction of the initial 
energy deposited in the initial object ("fireball") which is 
then the subject of further investigations. This fraction 
is called inelasticity and it is relevant also for low energy 
nuclear reactions [7]. 

The knowledge of the energy deposited in the central 
rapidity region in heavy ion collisions at RHIC and LHC 
is crucial [8] . Dividing this number by the volume of the 
formed system, we will have an estimate of the initial 
energy density in such collisions. If it is high enough we 
may be in a new phase of hadronic matter: the plasma 
of quarks and gluons (QGP). 

On the other hand, the knowledge of the momentum 
spectrum of the particles measured in the large rapid- 
ity region, and, in particular, those with the quantum 
numbers of the projectile (the so called leading parti- 
cles or LP) gives valuable information about the non- 
perturbative dynamics of QCD. Moreover, the LP spec- 
trum and the inelasticity of the reaction are very useful 
in cosmic ray physics, in the description of the evolution 
of hadronic showers in the atmosphere [9] . 

In the model considered here we hope to extract in- 
formation about the gluonic structure of hadrons from 
observablcs like mass, diffractive mass and leading par- 
ticle spectra, which are, at least in principle, very easy 
to measure. This model describes only certain aspects 
of hadronic collisions, related to energy flow and energy 
deposition in the central rapidity region. It should not be 
regarded as an alternative to a field-theoretical approach 
to amplitudes, but rather as an extension of the naive 
parton model. The reason for using it is that it may be 
good enough to account for energy flow in an economic 
way. The deeper or more subtle aspects of the underlying 



field theory probably (this is our belief) do not manifest 
themselves in energy flow, but rather in other quantities 
like the total cross section. Inspite of its simplicity, this 
model can teach us a few things and predict another few. 
This is encouraging because in the near future new data 
from FERMILAB, RHIC and LHC will be available. 

D. A brief history of the IGM 

Long time ago, based on qualitative ideas advanced 
by Pokorski and Van Hove [10], we started to develop 
a model to study energy deposition, connecting it with 
the apparent dominance of multiparticle production pro- 
cesses by the gluonic content of the impinging hadrons, 
hence its name: Interacting Gluon Model (IGM) [11]. Its 
original application to the description of inelasticity [12] 
and multiparticle production processes in hydrodynam- 
ical treatments [13] was followed by more refined appli- 
cations to leading charm production [14] and to single 
diffraction dissociation, both in hadronic reactions [15] 
and in reactions originated by photons [16]. These works 
allowed for providing the systematic description of the 
leading particle spectra [17] and clearly demonstrated 
that they are very sensitive to the amount of gluonic 
component in the diffracted hadron as observed in [18] 
and [19]. We have found it remarkable that all the results 
above were obtained using the same set of basic param- 
eters with differences arising essentially only because of 
the different kinematical limits present in each particular 
application. All this points towards a kind of universality 
of energy flow patterns in all the above mentioned reac- 
tions. The IGM was further developed and fluctuations 
in impact parameter were included in [20] , where a care- 
ful study of the inelasticity in proton-nucleus reactions 
was performed. The model was employed by the Camp- 
inas group of cosmic ray physics to reanalyse data from 
the AKENO collaboration and extract the proton-proton 
and proton-air cross sections [21]. This group used the 
IGM also to study the nucleonic and hadronic fluxes in 
the atmosphere [22]. 

Recent experimental developments encouraged us to 
return to the IGM picture of energy flow. One of them 
was connected with the new, more refined data on the 
leading proton spectra in ep — > e'pX obtained recently 
by the ZEUS collaboration [23] . Another one was a recent 
work on central mass production in Double Pomeron Ex- 
change (DPE) process reported in [24] allowing, in prin- 
ciple, for the extraction of the Pomeron-Pomeron total 
cross section oipip (see [25]). Finally, in the last years 
it became possible to study low x physics experimen- 
tally, at HERA and at RHIC. In this regime we probe 
the very low momentum region of the gluon distributions 
in hadrons and nuclei, where a qualitatively new behav- 
ior is expected to be dominant. This newly explored 
sector of the hadronic wave function is called by some 
authors Color Glass Condensate [26]. Its gluon density 
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is so high that the gluonic system can be treated semi- 
classically and in the weak coupling regime. This is the 
Bose-Einstein condensate of the strong interactions at 
the fundamental level and has been object of experimen- 
tal searches. Some of its signatures are related to energy 
flow observables and, in particular, to the leading parti- 
cle spectrum. Indeed, in [27] it was suggested that the 
LP spectrum in the saturation regime will go through a 
dramatically softening. In some earlier works we have 
predicted a slow softening but we did not include satu- 
ration. In view of the relevance of this subject, we plan 
to address this problem in the near future. 

In the next section we shall provide a brief description 
of the IGM, stressing the universality of energy flow and 
then we devote the other sections to discuss the applica- 
tions of the model. 



II. THE MODEL 

The IGM is based on the idea that since about half 
of a hadron momentum is carried by gluons and since 
gluons interact more strongly than quarks, during a high 
energy hadron-hadron collision there is a separation of 
constituents. Valence quarks tend to be fast forming 
leading particles whereas gluons tend to be stopped in 
the central rapidity region. The collision between the two 
gluonic clouds is treated as an incoherent sum of multi- 
ple gluon-gluon collisions, the valence quarks playing a 
secondary role in particle production. While this idea 
is well accepted for large momentum transfer between 
the colliding partons, being on the basis of some mod- 
els of minijet and jet production [1,28-34], in the IGM 
(and also in [29] and [34]) its validity is extended down to 
low momentum transfers, only slightly larger than Aqcd- 
At first sight this is not justified because at lower scales 
there are no independent gluons, but rather a highly cor- 
related configuration of color fields. There are, however, 
some indications coming from lattice QCD calculations, 
that these soft gluon modes are not so strongly corre- 
lated. One of them is the result obtained in [35] , namely 
that the typical correlation length of the soft gluon fields 
is close to 0.3 fm. Since this length is still much smaller 
than the typical hadron size, the gluon fields can, in a 
first approximation, be treated as uncorrelated. Another 
independent result concerns the determination of the typ- 
ical instanton size in the QCD vacuum, which turns out 
to be of the order of 0.3 fm [36]. As it is well known (and 
has been recently applied to high energy nuclcon-nucleon 
and nucleus- nucleus collisions) instantons are very im- 
portant as mediators of soft gluon interactions [37]. The 
small size of the active instantons lead to short distance 
interactions between soft gluons, which can be treated as 
independent. 

These two results taken together give support to the 
idea that a collision between two gluon clouds may be 
viewed as a sum of independent binary gluon-gluon colli- 



sions, which is the basic assumption of our model. De- 
veloping the picture above with standard techniques and 
enforcing energy-momentum conservation, the IGM be- 
comes the ideal tool to study energy flow in high energy 
hadronic collisions. Confronting this simple model with 
several and different data sets we obtained a surprisingly 
good agreement with experiment. 




FIG. 1. Schematic IGM pictures for (a) non-diffractive 
(ND), (b) and (c) single diffractive (SD) and (d) double 
Pomeron exchange (DPE) processes. 
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FIG. 2. Phase space limits of ND, SD and DPE processes 
in the IGM. The ^ In ^ line in a) indicates the rapidity Y of 
the produced mass M. 
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The hadron-hadron interaction follows the simple pic- 
ture shown in Fig. 1: the valence quarks fly through es- 
sentially undisturbed whereas the gluonic clouds of both 
projectiles interact strongly with each other forming a 
central fireball (CF) of mass M. The two incoming pro- 
jectiles pi and P2 loose fractions x and y of their original 
momenta and get excited forming what we call leading 
jets (LJ's) carrying xl = 1 — x and = 1 — y fractions 
of the initial momenta. Depending on the type of the 
process under consideration one encounters the different 
situations depicted in Fig. 1. 

In a non-diffractive (ND) process (Fig. la) one central 
fireball of mass M is formed, whereas in single diffractive 
(SD) events (Figs, lb and c) the corresponding diffrac- 
tive systems have masses Mx or My (comprising also the 
mass of CF) . In double Pomeron exchanges (DPE) (Fig. 
Id) a CF of mass Mxy is formed. In Fig. 2 we show 
their corresponding phase space limits. The only differ- 
ence between ND and SD or DPE processes is that in the 
latter cases the energy deposition is done by a restricted 
subset of gluons which in our language is a "kinemati- 
cal" Pomeron (IP) , the name which we shall use in what 
follows. 

The central quantity in the IGM is x( x j y)i the proba- 
bility to form a CF carrying momentum fractions x and y 
of two colliding hadrons. It follows from the quantitative 
implementation of the ideas described above. The essen- 
tial ingredients are the assumption of multiple indepen- 
dent gluon-gluon collisions, low momentum dominance 
of the gluon distributions and energy-momentum conser- 
vation. The derivation of our main formula, presented 
below, can be found in Appendix A. x(x, y) is given by: 

X(x,y) = — %= exp{--^[(2/ 2 )^- 

iTtyJ L) X y £L> X y 

+ (x 2 )(y-(y)) 2 +2(xy)(x-(x))(y-(y))}} (1) 
where D xy = (x 2 )(y 2 ) — (xy) 2 and 

, x ,n r^dy'ylm^X^y'), (2) 

Jo 

with xo defined by the normalization condition 



(x n y m ) = / dx' 
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[ dx [ dy x(x, y) 6{xy - K mm 
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1, (3) 



where K min — ^ is the minimal inelasticity defined by 
the mass mo of the lightest possible CF. The function 
ui(x',y'), sometimes called the spectral function, repre- 
sents the average number of gluon-gluon collisions as a 
function of x' e y'\ 



uj(x',y') 



dn 



dx' dy' 



(4) 



The appearance of the number n comes from the use 
of Poissonian distributions, which, in turn, is a conse- 
quence of the assumption of independent gluon-gluon 



collisions. u)(x',y') contains all the dynamical inputs of 
the model both in the perturbative (semihard) and non- 
perturbative (soft) regimes. The soft and semihard com- 
ponents are given by: 



1 X ' y >~ a ( S ) 

x 8(x'y'-K min )6(^^-x'y') (o) 



G(x')G(y') 
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and 



^{x',y') = 



a{s) 



G{x')G{y>) 



x 6{x'y' - 4PT """ ; 



(6) 



where a gg and a^ g are the soft and semihard gluonic cross 
sections, pr min is the minimum transverse momentum for 
minijet production and a denotes the impinging projec- 
tiles cross section. 

The values of x max and y m ax depend on the type of 
the process under consideration. For non-diffractive pro- 
cesses all phase space contained in the shaded area is 
allowed and in this case we have: 



Xmax — ymax 



1 



(7) 



The effective number of gluons from the correspond- 
ing projectiles are denoted by G's and have been approx- 
imated in all our works by the respective gluon distri- 
bution functions. There has been a remarkable progress 
in the knowledge of the parton distributions in hadrons 
[38-40], especially in the low x region, which becomes 
crucial at energies in the TeV range. Since in our pre- 
vious applications of the IGM we have been studying 
collisions in the GeV domain, there was no need to use 
very sophisticated parton distributions. Moreover, very 
often we needed parton densities at very low scales, which 
were not considered in the analyses presented in [38-40]. 
In some cases, we have used the parametrization of [41], 
which is better suited for small scales. However, as it will 
be shown, the IGM can describe both the hadronic and 
nuclear collision data with the following simple form of 
the gluon distribution function in the nucleon: 



G(x) = p(m + 1) 



(1 



(8) 



with m = 5 and the fraction of the energy-momentum 
allocated to gluons is equal to p — 0.5. 

In the IGM picture, diffractive and non-diffractive 
events have been treated on the same footing in terms 
of gluon-gluon collisions. Single diffractive processes re- 
ceive great attention mainly because of their potential 
ability to provide information about the most important 
object in the Regge theory, namely the Pomeron (IP), its 
quark-gluon structure and cross sections. As can be seen 
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in Fig. lb (lc), the "diffractivc mass" M x (M Y ) is just 
the invariant mass of a system composed of the CF and 
LJ formed by one of the colliding projectiles. The main 
difference with the "non-diffractive mass" M in Fig. la 
is that the energy transfer from the diffracted projectile 
is now done by the highly correlated subset of gluons (de- 
noted by IP) which are supposed to be in a color singlet 
state. In technical terms it means that in comparison to 
the previous applications of the IGM cited before, we are 
free to change both the possible shape of the function 
G(x) (= Gip(x)), the number of gluons participating in 
the process and the cross section a (= a v ip) in the spec- 
tral function oj used above. The function Gip(x) should 
not be confused with the momentum distribution of the 
gluons inside the Pomeron, f g /ip((3). The former is given 
by the convolution of the latter with the Pomeron flux 
factor as discussed in the Appendix B. Actually we have 
found that we can keep the shape of G(x) the same as 
before and the only change necessary to reach agreement 
with data is the amount of energy-momentum p = pjp 
allocated to the impinging hadron and which will find 
its way in the object that we call IP. It turns out that 
Pip ~ 0.05, whereas p ~ 0.5 for all gluons encountered 
so far. This choice, with m = 5 in cq. (8), corresponds 
to an intermediate between "soft" and "hard" Pomeron 
(see Appendix B) and will be used in what follows. Just 
in order to make use of the present knowledge about the 
Pomeron, we have chosen a(s) = a p ip = a + b ln(s/so) 
where So = 1 GeV 2 and a — 2.6 mb and b — 0.01 mb. 

In single diffractive processes only a limited part of the 
phase space supporting the x( x j y) distribution is allowed 
and in this case the integration limits in the moments of 
the spectral function w (eq. (2)) depend on the mass Mx 
or My that is produced: 

Xmax — 1 ! Umax V ] X max Umax — Mx/s (9) 
Xmax — X , Umax — 1 ; X max Umax — M-y / S (10) 

By reducing these maximal values we select events in 
which the energy released by the projectile emitting IP 
is small and at the same time allow the formation of a ra- 
pidity gap betwen the diffractive mass and the diffracted 
projectile. This is the experimental requirement defining 
a SD event. 

Double Pomeron Exchange processes, inspite of their 
small cross sections, are inclusive measurements and do 
not involve particle identification, dealing only with en- 
ergy flow. Such a process was recently measured by UA8 
[24] and used to deduce the IP IP cross section, app/p . It 
turned out that using this method one gets (Jipip which 
apparently depends on the produced mass Mxy ■ This 
fact was tentatively interpreted as signal of glueball for- 
mation [24]. In the IGM a double Pomeron exchange 
event (Fig. Id) is seen as a specific type of energy flow. 
The difference between it and the "normal" energy flow 
as represented by Fig. la is that now the gluons involved 
in this process must be confined to the object we called IP 



above. We are implicitly assuming that all gluons from 
Pi and pi participating in the collision (i.e., those emit- 
ted from the upper and lower vertex in Fig. Id) have to 
form a color singlet. In this case two large rapidity gaps 
will form separating the diffracted hadron p\ , the Mxy 
system and the diffracted hadron p2, which is the exper- 
imental requirement defining a DPE event. Also in this 
case only a limited part of the phase space supporting 
the x( x iV) distribution is allowed and the limits in the 
moments of the spectral function uj (eq. (2)) depend on 
the mass Mxy in the following way: 

Xmax — X J ymax — U 5 X max ]Jraax — M XY I ' S ' (H) 

As before Gip(x) represent the number of gluons par- 
ticipating in the process and the cross section a, appear- 
ing in eqs. (5) and (6), represents now the Pomeron- 
Pomeron cross section, opp/p. 

The clear separation between valence quarks and 
bosonic degrees of freedom does not appear exclusively in 
the IGM. It appears also in soliton models of the nucleon 
[42]. In the Chiral Quark Soliton Model [43], for exam- 
ple, the nucleon is made of three massive quarks bound 
by the self-consistent pion field (the "soliton"). It is in- 
teresting to observe that, according to this model, in a 
collision of two nucleons the valence quarks would inter- 
act much less than the pions and therefore would filter 
through and populate the large rapidity regions leaving 
behind a blob of pionic matter in the central region. 

III. INELASTICITY 

The energy dependence of inelasticity is an important 
problem which is still subject of debate [44]. Generally 
speaking, inelasticity K is the fraction of the total en- 
ergy carried by the produced particles in a given colli- 
sion. However in the literature one finds several possible 
ways to define it. In the first one, inelasticity is defined 
as 

Jt, = ^ (12) 

where y/s is the total reaction energy in its center of mass 
frame and M is the mass of the system (fireball, string, 
etc.) which decays into the final produced particles. The 
second definition of K considered here is 

k 2 = J dy ^~^j coshy ( 13 ) 



where fii = \jp% i + mf is the transverse mass of the pro- 
duced particles of type i and drii/dy their measured ra- 
pidity distribution. These two definitions are, in princi- 
ple, model independent, although the mass M might be 
difficult to evaluate in certain models. 

The main difference between K\ and K 2 is that, 
whereas the first one refers to partons, the second one 



5 



refers to final observed hadrons. K2 implicitly includes 
the kinetic energy of the object of mass M. From the 
theoretical point of view, K\ is a very interesting quan- 
tity because it can be easy to calculate and because it 
is the relevant quantity when studying the formation of 
dense systems (e.g. quark-gluon plasma). 

In Ref. [12] we used the IGM to study the energy de- 
pendence of K\ . We concluded that the introduction of 
a semihard component (minijets) in that model produces 
increasing inelasticities at the partonic level. In Ref. [13] 
we introduced a hadronization mechanism in the IGM, 
calculated the rapidity distributions of the produced par- 
ticles, compared our results with the UA5 and UA7 data 
and finally calculated K-i- The purpose of this exercise 
was to verify whether the hadronization process changed 
our previous conclusion. We found that, whereas some 
quantitative aspects, like the existence or not of Feyn- 
man Scaling in the fragmentation region and the numer- 
ical values of K 2 , depend very strongly on details of the 
fragmentation process, the statement that minijets lead 
to increasing inelasticities remains valid. 

In Fig. 3 we show the IGM pseudorapidity distribu- 
tions compared to UA5 data at different energies [45] 
and CDF [46] data at y/s = 1800 GeV. For the sake of 
comparison with other models based both on soft and 
semihard dynamics, we show in Fig. 4 our results for the 
multiplicity (Fig. 4a) and central rapidity density (Fig. 
4b) together with the results of HIJING [1] for the same 
quantities. 




FIG. 3. Pseudorapidity distributions measured at the cen- 
tral rapidity region. Data are from the UA5 collaboration 
[45] at different energies and from CDF collaboration [46] at 
<fs = 1800 GeV. Full lines show the IGM results. 



Both models fit the data but differ significantly when 
one switches off the semihard (minijet) contribution. 
Whereas in HIJING Feynman Scaling violation in the 
central region (the growth of dn/dr)\ with y/s) is 



entirely due to the minijets, in the IGM this behavior is 
partly due to soft interactions, there being only a quan- 
titative difference when minijets are included. 

In Fig. 5 we plot K2 (full lines) and K\ (dashed lines 
as a function of yfs. The lower curves show the results 
when minijets are switched off and only soft interactions 
take place. The upper curves show the effect of including 
minijets. 




ol 1 1 1 1 

10 z 10 3 10" 




FIG. 4. a) Average charged multiplicities as a function of 
the reaction energy. Squares and circles are experimental 
data. Full lines show the IGM results with and without the 
semihard contribution (lower curve). Dashed lines show the 
same quantities calculated with HIJING [1]. b) The same as 
a) for the central pseudorapidity distribution dn/drj\ . 



IV. LEADING CHARM AND BEAUTY 

In Ref. [14] we treated leading charm production in 
connection with energy deposition in the central rapid- 
ity region giving special attention to the correlation be- 
tween production in central and fragmentation regions. 
The significant difference between the xp dependence of 
leading and nonleading charmed mesons [47] was not pos- 
sible to be explained with the usual perturbative QCD 



G 



[48] or with the string fragmentation model contained in 
PYTHIA [49]. 
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FIG. 5. Inelasticities K\ (dashed lines) and K2 (solid lines) 
with minijets (upper curves) and without minijets (lower 
curves) as a function of the reaction energy. 



In the case of pion-nucleon scattering, the measured 
leading charmed mesons [47] are D~ and the nonleading 
are D + . In the spirit of IGM, the central production ig- 
nores the valence quarks of target and projectile which, 
in the first approximation, just "fly through" . Because of 
this, the centrally produced D's do not show any leading 
particle effect. There are, however, two distinct ways to 
produce D mesons out of LJ's: fragmentation and recom- 
bination. We assumed that, whenever energy allows, we 
would have also cc pairs in the LJ (produced, for example, 
from the remnant gluons present there). These charmed 
quarks might undergo fragmentation into D mesons and 
also recombine with the valence quarks. Whereas D + 
and D~ mesons are equally produced via fragmentation 
only "leading" D~'s (which carry the valence quarks of 
target and projectile) can be produced by recombination. 
It turns out that only this last process will produce asym- 
metry. 

The idea that cc pairs pre-exist in the projectiles can 
be made more precise and the origin of these "intrinsic 
charm" pairs can be attributed to the existence of a me- 
son cloud around the nucleons and pions [50] . 

The asymmetry in D meson production can be defined 
as: 



A{x F ) 



da D (x F ) 
dx p 



d(T D ^ (xp) 



"(*F) 



(14) 



dXF 



dXF 



In Fig. 6 we compare our calculations with experi- 
mental data from the WA82, E769 and E791 [47] collab- 
orations. The main conclusion of the work is that if one 
takes properly into account the correlation between energy 
deposition in the central region and the leading particle 
momentum distribution, at higher energies the increase 
of inelasticity will lead to the decrease of the asymmetry 



in heavy quark production. In other words, if the fraction 
of the reaction energy released in the central region in- 
creases the asymmetry in the x p distributions of charmed 
mesons will become smaller. In Fig. 7 we illustrate this 
quantitatively and also consider the leading beauty pro- 
duction. 




s = 26GeV 



FIG. 6. Asymmetry calculated with the IGM and com- 
pared with WA82 (solid circles), with E769 (open squares) 
and E791 (open triangles) data [47]. Solid, dashed and dot- 
ted lines correspond to different weights of the recombination 
component equal to 80%, 50% and 20% respectively. 
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FIG. 7. B~/B + (solid lines) and D~/D+ (dashed lines) 
asymmetries at ^/s = 26 and 1800 GeV. 



V. DIFFRACTIVE MASS SPECTRA IN 
HADRON-HADRON COLLISIONS 

Diffractive scattering processes have received increas- 
ing attention for several reasons. They are also related 
to the large rapidity gap physics and the structure of 
the Pomeron. In a diffractive scattering, one of the in- 
coming hadrons emerges from the collision only slightly 
deflected and there is a large rapidity gap between it and 
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the other final state particles resulted from the other ex- 
cited hadron. Diffraction is due to the Pomeron exchange 
but the exact nature of the Pomeron in QCD is not yet 
elucidated. The first test of a model of diffractive disso- 
ciation (SD) is the ability to properly describe the mass 
(Mx) distribution of diffractive systems, which has been 
measured in many experiments [51] and parametrized as 
(M x )~ a with a ~ 1. 

In Ref. [15] we studied diffractive mass distributions 
using the Interacting Gluon Model focusing on their en- 
ergy dependence and their connection with inelasticity 
distributions. One advantage of the IGM is that it was 
designed in such a way that the energy-momentum con- 
servation is taken care of before all other dynamical as- 
pects. This feature makes it very appropriate for the 
study of energy flow in high energy hadronic and nuclear 
reactions. As mentioned before, in our approach the def- 
inition of the object IP (see Fig. lb and c) is essentially 
kinematical, very much in the spirit of those used in other 
works which deal with diffractive processes in the parton 
and/or string language. In order to regard our process 
as being of the SD type we simply assume that all gluons 
from the target hadron participating in the collision (i.e., 
those emitted from the lower vertex in Fig. lb) have to 
form a colour singlet. Only then a large rapidity gap 
will form separating the diffracted hadron and the Mx 
system. Otherwise a colour string would develop, con- 
necting the diffracted hadron and the diffractive cluster, 
and would eventually decay, filling the rapidity gap with 
produced secondaries. As was said above, it is a special 
class of events in which the energy released by the projec- 
tile emiting IP is small and consequently the diffractive 
mass is small. Once only a limited part of the phase 
space is allowed, the limits in the moments of the spec- 
tral function ui (eq. (2)), depend on the mass Mx that is 
produced through the constraint y max = y = M x / s (see 
eq. (9)). 
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FIG. 8. Diffractive mass spectrum for pp collisions cal- 
culted with the IGM and compared with CERN-ISR data 
[52]. 



As shown in Appendix A, the mass spectra for SD 
processes is given by: 



1 da 



adM x 



dN 
dM x 



dx 





sy) 



dyx(x,y) 



(xy - K„ 



(15) 



Although in the final numerical calculations the above 
complete formulation is used, it is worthwhile to present 
approximate analytical results in order to illustrate the 
main characteristic features of the IGM diffractive dis- 
sociation processes. By keeping only the most singular 
terms in gluon distribution functions, i.e., G(x) ~ ^ and 
only the leading terms in y/s, as shown in Appendix A, 
we arrive at the following expression: 



dN 
dM\ 



±H(M x ,s)F(M x ,s) 



const 
~M\ 



exp 



c In 



c In 



Mi. 



In 



Mi 



(16) 



where c is a constant and mo is a soft energy scale. The 
expression above is governed by the -7^5- term. The other 

two terms have a weaker dependence on M x . They dis- 
tort the main ( -^jj- ) curve in opposite directions and tend 

to compensate each other. It is therefore very interesting 
to note that even before choosing a very detailed form 
for the gluon distributions and hadronic cross sections 
we obtain analytically the typical shape of a diffractive 

spectrum, -ttt. 

x 

In Fig. 8 we show our diffractive mass spectrum and 
compare it to experimental data from the CERN-ISR 
[52] . Fig. 9 shows the diffractive mass spectrum for yfs — 
1800 GeV compared to experimental data from the E710 
Collaboration [53]. In these curves we have used our 
intermediate Pomeron profile: Gip(y) given by (8) with 
ttl — 5 and pjp = 0.05. 

In all curves we observe a modest narrowing as the en- 
ergy increases. This small effect means that the diffrac- 
tive mass becomes a smaller fraction of the available en- 
ergy t/s. In other words, the "diffractive inelasticity" de- 
creases with energy and consequently the "diffracted lead- 
ing particles" follow a harder xp spectrum. Physically, 
in the context of the IGM, this means that the deposited 
energy is increasing with y/s but it will be mostly released 
outside the phase space region that we are selecting. A 
measure of the "diffractive inelasticity" is the quantity 
£ = M x /s. It is very simple to calculate its average 
value (£) from the diffractive mass spectrum. Making a 
trivial change of variables we get: 



(0 (a) 



- dN 



(17) 



where £ m i„ (= 1.5/ s) and £ maa; (= 0.1) are the same used 
in other works. In Fig. 10 we plot (£) against ^fs. As 
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it can be seen (£) decreases with y/s not only because 
Cmin becomes smaller but also because dN/d£ changes 
with the energy, falling faster. Also shown in Fig. 10 is 
the quantity (sometimes used in connection with the 
energy dependence of the single diffractive cross-section) 
for e — 0.08 (dashed lines) and e = 0.112 (dotted lines). 
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FIG. 9. Diffractive mass spectrum for pp collisions cal- 
culted with the IGM and compared with FERMILAB Teva- 
tron data [53]. 
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VI. DIFFRACTIVE MASS SPECTRA IN 
ELECTRON-HADRON COLLISIONS 

At the HERA electron-proton collider the bulk of the 
cross section corresponds to photoproduction, in which 
a beam electron is scattered through a very small angle 
and a quasi-real photon interacts with the proton. For 
such small virtualities the dominant interaction mech- 
anism takes place via fluctuation of the photon into a 
hadronic state (vector meson dominance) which interacts 
with the proton via the strong force. High energy pho- 
toproduction therefore exhibits similar characteristics to 
hadron-hadron interactions. 

In Ref. [16] we studied diffractive mass distributions in 
a photon-proton collision. The photon is converted into 
a mesonic state and then interacts with the incoming 
proton. The diffractive meson-proton interaction follows 
then the usual IGM picture. The diffracted proton in 
Fig. lb), looses only a fraction y of its momentum but 
otherwise remains intact. In the limit y — > 1, the whole 
available energy is stored in Mj which then remains at 
rest, i.e., Yx = 0. For small values of y we have small 
masses Mx located at large rapidities Yx- As before 
the upper cut-off y m ax (= V = M x /s) is a kinemat- 
ical restriction preventing the gluons coming from the 
diffracted proton (and forming our object IP) to carry 
more energy than what is released in the diffractive sys- 
tem. It plays a central role in the adaptation of the IGM 
to diffractive dissociation processes being responsible for 
its proper M x dependence. 
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FIG. 11. Diffractive mass spectrum for 7p collisions at 
W = 187 GeV calculted with the IGM and compared with 
HI data [54]. The different curves correspond to the choices: 
I (m = 0.31 GeV, a = 2.7 mb), II (m = 0.35 GeV, 
a = 2.7 mb), III (m = 0.31 GeV, a = 5.4 mb) and IV 
(mo = 0.35 GeV, a — 5.4 mb), respectively. 



In the upper leg of Fig. lb) we have assumed, for sim- 
plicity, the vector meson to be p° and take G p (x) = 
G"{x). Since the parameter p/a appearing in eqs. (5) 
and (6) has been fixed considering the proton-proton 
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diffractive dissociation and here we are addressing the 
pp° case there exists some freedom to change a. We can 
also investigate the effect of small changes in the value of 
mo on our final results. 
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FIG. 12. Diffractive mass spectrum for 7p collisions at 
W = 187 GeV calculted with the IGM and compared with 
HI data [54]. The solid line (curve II) corresponds to the 
choice mo = 0.35 GeV, a — 2.7 mb and Gip(y). Curves I 
(dashed) and III (dotted) are obtained replacing Gip(y) by 
Gjp{y) and G B IP (y) respectively. Curve IV is obtained with 
G s IP (y) and mo = 0.50 GeV and a = 5.4 mb. 

In Fig. 11 we compare our results, eq. (15), for dif- 
ferent choices of too and a with the data from the HI 
collaboration [54]. In all these curves we have used our 
intermediate Pomeron profile. In Fig. 12 we compare the 
same data with our mass spectrum obtained with G l } P (y) 
(curve I), Gip(y) (curve II) and Gj P (y) (curve III). This 
comparison suggests that the "hard" Pomeron can give a 
good description of data. The same can be said about our 
"mixed" Pomeron, which, in fact seems to be more hard 
than soft. These three curves were calculated with ex- 
actly the same parameters and normalizations, the only 
difference being the Pomeron profile. Soft and hard gluon 
distributions Gj P (y) are calculated in the Appendix B. 
Apparently the "soft" Pomeron (curve III) is ruled out 
by data. 



VII. DIFFRACTIVE MASS SPECTRA IN 
DOUBLE POMERON EXCHANGE 

After ten years of work at HERA, an impressive 
amount of knowledge about the Pomeron has been accu- 
mulated, especially about its partonic composition and 
parton distribution functions. Less known are its inter- 
action properties. Whereas the Pomeron- nucleon cross 
section has been often discussed in the literature, the 
recently published data by the UA8 Collaboration [24] 
have shed some light on the Pomeron-Pomeron interac- 
tion. In [24] the Double Pomeron Exchange cross sec- 
tion was written as the product of two flux factors with 
the IP IP cross section, o~ipip, being thus directly pro- 
portional to this quantity. This simple formula relies on 



the validity of the Triple-Regge model, on the universal- 
ity of the Pomeron flux factor and on the existence of 
a factorization formula for DPE processes. However, for 
these processes the factorization hypothesis has not been 
proven and is still matter of debate [55-58] . In [59] it was 
shown that factorizing and non-factorizing DPE models 
may be experimentally distinguished in the case of dijet 
production. 

Fitting the measured mass spectra allowed for the de- 
termination of a ip ip and its dependence on Mx, the 
mass of the diffractive system. The first observation 
of the UA8 analysis was that the measured diffractive 
mass (Mx) spectra show an excess at low values that 
can hardly be explained with a constant (i.e., indepen- 
dent of Mx) (tip ip- Even after introducing some mass 
dependence in o~ipip they were not able to fit the spec- 
tra in a satisfactory way. Their conclusion was that the 
low Mx excess may have some physical origin like, for 
example glueball formation. 

Although the analysis performed in [24] is standard, 
it is nevertheless useful to confront it with the IGM de- 
scription of the diffractive interaction. Double Pomeron 
exchange processes, inspite of their small cross sections, 
appear to be an excellent testing ground for the IGM be- 
cause they are inclusive measurements and do not involve 
particle identification, dealing only with energy flow. In 
Ref. [25] we studied the diffractive mass distribution ob- 
served by UA8 Collaboration in the inclusive reaction 
pp -> pXp at y/s = 630 GeV, using the IGM with DPE 
included. The interaction follows the picture shown in 
Fig. Id. 

As shown in Appendix A, the mass spectrum for DPE 
processes is given by: 



1 da 



dN 



a dM X Y dM X Y 
x S(M X y 




dyx(x,y) 



K 2 



(18) 



As indicated in the recent literature [55-59] , one of the 
crucial issues in diffractive physics is the possible break- 
down of factorization. As stated in [56] one may have 
Regge and hard factorization. Our model does not rely 
on any of them. In the language used in [56], we need 
and use a "diffractive parton distribution" and we do not 
really need to talk about "flux factor" or "distribution of 
partons in the Pomeron" . Therefore there is no Regge 
factorization implied. However, we will do this connec- 
tion in eq. (B6) of the Appendix B, in order to make 
contact with the Pomeron pdf 's parametrized by the HI 
and ZEUS collaborations. As for hard factorization, it is 
valid as long as the scale p is large. In the IGM, as it will 
be seen, the scale is given by p 2 — xys, a number which 
sometimes is larger than 3 — 4 GeV 2 but sometimes is 
smaller, going down to values only slightly above Aq CD . 
When the scale is large (p 2 > p% . ) we employ Eq. (6) 
and when it is smaller (to,q < p? < p\, ) we use Eq. 
(5). Therefore, in part of the phase space we are inside 
the validity domain of hard factorization, but very often 
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we are outside this domain. From the practical point of 
view, Eq. (6), being defined at a semihard scale, relies 
on hard factorization for the elementary gg — > gg inter- 
action, uses parton distribution function extracted from 
DIS and an elementary cross section a gg taken from stan- 
dard pQCD calculations. The validity of the factorizing- 
like formula Eq. (5) is an assumption of the model. In 
fact, the relevant scale there is m$ ~ ^qcd an< ^' strictly 
speaking, there are no rigorously defined parton distribu- 
tions, neither elementary cross sections. However, using 
Eq. (5) has non-trivial consequences which were in the 
past years supported by an extensive comparison with 
experimental data. In this approach, since we have fixed 
all parameters using previous data on leading particle 
formation and single diffractive mass spectra, there are 
no free parameters here, except uipip ■ 

We start evaluating Eq. (18) with the inputs that were 
already fixed by other applications of the IGM, namely, 
(B7) with pd — 0.05. In Fig. 13 we show the numeri- 
cal results for DPE mass distribution normalized to the 
"AND" data sample of [24] . We have fixed the parameter 
a (= o-pip) appearing in Eq. (5) and Eq. (6), to 0.5mb 
(dashed lines) and 1.0 mb (solid lines). 

We emphasize that, in this approach, since we have 
fixed all parameters using previous data on leading par- 
ticle formation and single diffractive mass spectra, there 
are no free parameters here, except app. As it can be 
seen from the figure, in our model we obtain the fast in- 
crease of spectra in the low mass region without the use 
of a Mx dependent IP IP cross section and this quantity 
seems to be approximately o~pp ~ 0.5 mb. 
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FIG. 13. IGM DPE diffractive mass distributions: Solid 
and dashed lines show the results with aipip equal to 0.5 mb 
and 1.0 mb, respectively, calculated with the intermediate 
Pomeron profile. Our curves were normalized to the "AND" 
data sample of [24]. 
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FIG. 14. IGM DPE diffractive mass distributions: Solid 
line as in Fig. 13, dashed and dash-dotted lines represent 
the "hard" and "super-hard" Pomeron profiles. In all cases 
oipip — 0.5 mb. Our curves were normalized to the "AND" 
data sample of [24]. 

We next replace (B7) by the convolution (B6) to see 
which of the previously considered Pomeron profiles, hard 
or superhard, gives the best fit of the UA8 data. In doing 
so, we shall keep everything else the same, i.e., pd = 0.05 
and o~pp = 0.5 mb. 

In Fig. 14, we repeat the fitting procedure used in Fig. 
13 for these Pomeron profiles. Solid, dashed and dash- 
dotted lines represent respectively Eq. (B7), hard and 
superhard Pomerons. We see that, for harder Pomeron 
profiles we "dig a hole" in the low mass region of the 
spectrum. Note that the solid lines are the same as in 
Fig. 13. Looking at the figure, at first sight, we might 
be tempted to say that Eq. (B7) gives the best agree- 
ment with data and a somewhat worse description can 
be obtained with the hard Pomeron (in dashed lines), 
the superhard being discarded. However, comparing the 
dashed lines in Fig. 13 and Fig. 14 and observing that 
they practically coincide with each other, we conclude 
that the same curve can be obtained either with (B7) 
and <Jpp = 1.0 mb (dashed line in Fig. 13) or with (B6), 
(B2) and app = 0.5 mb (dashed line in Fig. 14). In other 
words we can trade the "hardness" of the Pomeron with 
its interaction cross section. The following two objects 
give an equally good description of data: i) a Pomeron 
composed by more and softer gluons and with a larger 
cross section and ii) a Pomeron made by fewer, harder 
gluons with a smaller interaction cross section. We have 
checked that this reasoning can be extended to the su- 
perhard Pomeron. Although, apparently disfavoured by 
Fig. 14 (dash-dotted lines), it might still fit the data 
provided that app < 0.25 mb. Given the uncertainties 
in the data and the limitations of the model, we will not 
try for the moment to refine this analysis. It seems pos- 
sible to describe data in a number of different ways. We 
conclude then that nothing exotic has been observed and 
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also that the Pomeron-Pomeron cross section is bounded 
to be app < 1.0 mb. 



10 



da , 



dM, 



(mb/ GeV) 



10 • 



10 ' 



10 




200 



50 100 150 

M x {GeV) 

FIG. 15. IGM prediction for da/dMx at LHC with 
a ip ip = 1.0 mb. Cross(+) and Cross(x) are predictions made 
by Brandt et al. [24] for two values of effective Pomeron in- 
tercepts a(0) = 1 + e. 

In Fig. 15 we compare our predictions for da/dMx 
(mb/GeV) for LHC (^/s = 14 TeV) assuming an M x - 
independent app — 1.0 mb (and using (B7)) with pre- 
dictions made by Brandt et al. [24] for two values of ef- 
fective Pomeron intercepts (a(0) = 1 + e), e = 0.0 and 
0.035. 
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FIG. 16. Ratio double/single Pomeron exchange mass dis- 
tributions as a function of Mx- In both cases we have as- 
sumed app — 1.0 mb (for DPE processes) and a v p = 1.0 mb 
(for SPE processes). 



Although the normalization of our curves is arbitrary, 
the comparison of the shapes reveals a difference be- 
tween the two predictions. Whereas the points (from 
[24]) show spectra broadening with the c.m.s. energy, we 
predict (solid line) the opposite behavior: as the energy 
increases we observe a (modest) narrowing for da/dMx- 
This small effect means that the diffractive mass becomes 
a smaller fraction of the available energy s/s. In other 
words, the "double diffractive inelasticity" decreases with 
energy in the same way as the "diffractive inelasticity" , 
as seen in Fig. 10. 

We are not able to make precise statements about the 
diffractive cross section (in particular about its normal- 
ization) with our simple model. Nevertheless, the nar- 
rowing of da d p e / dMx suggests a slower increase (with 
yfs) of the integrated distribution a pipe- We found this 
same effect also for asPE- This trend is welcome and is 
one of the possible mechanisms responsible for the sup- 
pression of diffractive cross sections at higher energies 
relative to some Regge theory predictions. 

In Fig. 16 we show the ratio R(Mx) defined by: 



R{M X ) 



1 da dpe 
&dpe dMx 

1 dasPE 
c spe dM x 



(19) 



This quantity involves only distributions previously nor- 
malized to unity and does not directly compare the cross 
sections (which are numerically very different for DPE 
and single diffraction). In R the dominant \/M\ fac- 
tors cancel and we can better analyse the details of the 
distributions which may contain interesting dynamical 
information. The most prominent feature of Fig. 16 is 
the rise of the ratio with Mx, almost by one order of 
magnitude in the mass range considered. This can be 
qualitatively attributed to the fact that, in single diffrac- 
tive events the object X has larger rapidities than the 
corresponding cluster formed in DPE events. As a con- 
sequence, when energy is released from the incoming par- 
ticles in a SPE event, it goes more to kinetic energy of 
the X system (i.e., larger momentum Px and rapidity 
Yx) and less to its mass. In DPE, although less energy 
is released, it goes predominantly to the mass Mx of the 
difractive cluster, which is then at lower values of Yx- 
In order to illustrate this behavior, we show in Fig. 17 
the rapidity distributions of the X (which has mass Mx) 
and XY (which has mass Mxy) systems. All curves 
are normalized to unity and with them we just want to 
draw attention to the dramatically different positions of 
the maxima of these distributions. The solid and dashed 
lines show 1/ada/dYx for DPE (curves on the left) and 
SPE (curves on the right) computed at y/s = 630 GeV 
and y/s — 2000 GeV, respectively. We can clearly observe 
that DPE and SPE rapidity distributions are separated 
by three units of rapidity and this difference stays nearly 
constant as the c.m.s. energy increases. The location of 
maxima in 1/ada/dYx and their energy dependence are 
predictions of our model. 
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FIG. 17. Double and single Pomeron exchange normalized 
rapidity (Yx) distributions. In both cases we have assumed 
app — 1.0 mb (for DPE processes) and a p p = 1.0 mb (for 
SPE processes). 



VIII. LEADING PARTICLE SPECTRA 

The leading particle effect is one of the most interest- 
ing features of multiparticle production in hadron-hadron 
collisions. In high energy hadron-hadron collisions the 
momentum spectra of outgoing particles which have the 
same quantum numbers as the incoming particles, also 
called leading particle (LP) spectra, have been measured 
already some time ago [60,61]. Later on, new data on 
pion-proton collisions were released by the EHS/NA22 
collaboration [62] in which the spectra of both outcoming 
leading particles, the pion and the proton, were simulta- 
neously measured. More recently data on leading pro- 
tons produced in eletron-proton reactions at HERA with 
a c.m.s. energy one order of magnitude higher than in 
the other above mentioned hadronic experiments became 
available [63] . In the case of photoproduction, data can 
be interpreted in terms of the Vector Dominance Model 
[64] and can therefore be considered as data on LP pro- 
duction in vector meson-proton collisions. These new 
measurements of LP spectra both in hadron-hadron and 
in eletron-proton collisions have renewed the interest on 
the subject, specially because the latter are measured at 
higher energies and therefore the energy dependence of 
the LP spectra can now be determined. 

It is important to have a very good understanding of 
these spectra for a number of reasons. They are the input 
for calculations of the LP spectra in hadron- nucleus colli- 
sions, which are a fundamental tool in the description of 
atmospheric cascades initiated by cosmic radiation [9,65]. 
There are several new projects in cosmic ray physics in- 
cluding the High Resolution Fly's Eye Project, the Tele- 
scope Array Project and the Pierre Auger Project [66] 
for which a precise knowledge of energy flow (LP spectra 
and inelasticity distributions) in very high energy colli- 
sions would be very useful. 



In a very different scenario, namely in high energy 
heavy ion collisions at RHIC, it is very important to know 
where the outgoing (leading) baryons are located in mo- 
mentum space. If the stopping is large they will stay in 
the central rapidity region and affect the dynamics there, 
generating, for example, a baryon rich equation of state. 
Alternatively, if they populate the fragmentation region, 
the central (and presumably hot and dense) region will 
be dominated by mesonic degrees of freedom. The com- 
position of the dense matter is therefore relevant for the 
study of quark gluon plasma formation [8] . 

In any case, before modelling p — A or A — A colli- 
sions one has to understand properly hadron-hadron pro- 
cesses. The LP spectra are also interesting for the study 
of diffractive reactions, which dominate the large xp re- 
gion. 

Since LP spectra are measured in reactions with low 
momentum transfer and go up to large xp values, it is 
clear that the processes in question occur in the non- 
perturbative domain of QCD. One needs then "QCD in- 
spired" models and the most popular are string models, 
like FRITIOF, VENUS or the Quark Gluon String Model 
(QGSM). Calculation of LP spectra involving these mod- 
els can be found in Refs. [67] and [68]. 

A. Leading particles in hadron-hadron collisions 

In the framework of the QCD parton model of high 
energy collisions, leading particles originate from the 
emerging fast partons of the collision debris. There is 
a large rapidity separation between fast partons and sea 
partons. Fast partons interact rarely with the surround- 
ing wee partons. The interaction between the hadron 
projectile and the target is primarily through wee par- 
ton clouds. A fast parton or a coherent configuration of 
fast partons may therefore filter through essentially unal- 
tered. Based on these observations and aiming to study 
p — A collisons, the authors of Ref. [67] proposed a mech- 
anism for LP production in which the LP spectrum is 
given by the convolution of the parton momentum dis- 
tribution in the projectile hadron with its correspond- 
ing fragmentation function into a final leading hadron. 
This independent fragmentation scheme is, however, not 
supported by leading charm production in pion-nucleus 
scattering. It fails specially in describing the D~/D + 
asymmetry. A number of models addressed these data 
and the conclusion was that valence quark recombination 
is needed. Translated to leading pion or proton produc- 
tion this means that what happens is rather a coalescence 
of valence quarks to form the LP and not an indepen- 
dent fragmentation of a quark or diquark to a pion or 
a nucleon. Another point is that the coherent config- 
uration formed by the valence quarks may go through 
the target but, due to the strong stopping of the gluon 
clouds, may be significantly decelerated. This correla- 
tion between central energy deposition due to gluons and 
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leading particle spectra was shown to be essential for the 
undertanding of leading charm production [14] . 

We follow the same general ideas of Ref. [67] but with a 
different implementation. In particular we replace inde- 
pendent fragmentation by valence quark recombination 
and free leading parton flow by deceleration due to "gluon 
stripping" . 

We have studied all measured LP spectra including 
those measured at HERA. We have found some universal 
aspects in the energy flow pattern of all these reactions. 
Universality means, in the context of the IGM, that 
the underlying dynamics is the same both in diffractive 
and non-diffractive LP production and both in hadron- 
hadron and photon-hadron processes. 

In Ref. [17] we analyzed leading particle spectra in 
hadronic collisions and, assuming VDM, the leading pro- 
ton spectra in e — p reactions. We have also considered 
the contribution coming from the diffractive processes. 
The leading particle can emerge from different regions 
of the phase space, according to the values assumed by 
%max and Vmax m eqs. (5) and (6). The distribution of 
its momentum fraction xl is given by: 

F{xl) = (1 - a) F nd (xL 

where a — ct\ 4- ct2 is the total fraction of single diffrac- 
tive (d) events from the lower and upper legs in Fig. 1, 
respectively. 
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FIG. 18. Comparison of our LP spectra F(xl) with data 
from [61] and [60]. 



Notice that a is essentially a new parameter here, 
which should be of the order of the ratio between the 
total diffractive and total inelastic cross sections [15]. 

In Fig. 18 we present our spectra of leading protons, 
pions and kaons respectively. The dashed lines show the 
contribution of non-diffractive LP production and the 
solid lines show the effect of adding a diffractive compo- 
nent, calculated with the intermediate Pomeron profile. 
All parameters were fixed previously and the only one 
to be fixed was a. For simplicity we have neglected the 
third diagram in Fig. 1 (c) , because it gives a curve which 
is very similar in shape to the non-diffractive curve. In 
contrast, the Pomeron emission by the projectile (Fig. 
lb) produces the diffractive peak. We have then chosen 
OL2 — and ol\ = a = 0.3 in all collision types. 

As expected, the inclusion of the diffractive component 
flattens considerably the final LP distribution bringing 
it to a good agreement with the available experimental 
data [60,61]. In our model there is some room for changes 
leading to fits with better quality. We could, for exam- 
ple, use a prescription for hadronization (as we did before 
in [13])) giving a more important role to it, as done in 
Ref. [67]. In doing this, however, we loose simplicity and 
the transparency of the physical picture, which are the 
advantages of the IGM. We prefer to keep simplicity and 
concentrate on the interpretation of our results. In first 
place it is interesting to observe the good agreement be- 
tween our curve and data for protons (Fig. 18a) in the 
low xl region. The observed protons could have been 
also centrally produced, i.e., they could come from the 
CF. However we fit data without the CF contribution. 
This suggests, as expected, that all the protons in this 
Xl range are leading, i.e., they come from valence quark 
recombination. In Figs. 18b) and 18c) we observe an 
excess at low xl- This is so because pions and kaons 
are light and they can more easily be created from the 
sea (centrally produced). Our distributions come only 
from the leading jet and consequently pass below the 
data points. A closer look into the three dashed lines 
in Fig. 18 shows that pion and kaon spectra are softer 
than the proton one. The former peak at x ~ 0.56 while 
the latter peaks at x ~ 0.62. In the IGM this can be 
understood as follows. The energy fraction that goes 
to the central fireball, K = Jxy^ is controled by the be- 
haviour of the function xi x i V) nd > which is approximately 
a double gaussian in the variables x and y, as it can be 
seen in expression (1). The quantities (x) and (y) play 
the role of central values of this gaussian. Consequently 
when (x) or (y) increases, this means that the energy 
deposition from the upper or lower leg (in Fig. 1) in- 
creases respectively. The quantities (x) and (y) are the 
moments of the uj function and are directly proportional 
to the gluon distribution functions in the projectile and 
target and inversely proportional to the target-projectile 
inelastic cross section. In the calculations, there are two 
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changes when we go from p — p to ir — p: 

(i) The first is that we replace crf^el by a inei wmcrl is 
smaller. This leads to an overall increase of the 
energy deposition. There are some indications that 
this is really the case and the inelasticity in tt — p 
is larger than in p — p collisions 1 . 

(n) The second and most important change is that 
we replace one gluon distribution in the proton 
G p {y) by the corresponding distribution in the pion 
G n (y). We know that G p (y) ~ (1 — y) 5 /y whereas 
G K {y) ~ (1 — y) 2 /y, i.e., that gluons in pions are 
harder than in protons. This introduces an asym- 
metry in the moments (a;) and (y) , making the lat- 
ter significantly larger. 

As a consequence, because of their harder gluon distri- 
butions, pions will be more stopped and will emerge from 
the collision with a softer xl spectrum. This can already 
be seen in the data points of Fig. 18. However since these 
points contain particles produced by other mechanisms, 
such as central and diffractive production, it is not yet 
possible to draw firm conclusions. One should mention 
here that there is another possible difference between nu- 
cleons and mesons which can contribute to the different 
behaviour of the leading particles in both cases. It is con- 
nected with the triple gluon junction present in baryons 
but not in mesons, which, if treated as an elementary 
object, can influence sunbstantially LP spectra (cf. [71]). 
We shall not discuss this possibility in this paper. 

The analysis of the moments (x) and (y) can also be 
done for the diffractive process shown in Fig. lb). Be- 
cause of the cuts in the integrations in eq. (2), they 
will depend on ix = 1 — y. We calculate them for 
p + p — ► p + X and 7r + p — ► 7r + X reactions. For low Xl 
they assume very similar values as in the non-diffractive 
case. For large xl however we find that (x) p ~ (x)^ and 
(y) p ~ (y)-K- The reason for these approximate equalities 
is that in diffractive processes we cut the large y' region 
and this is precisely where the pion and the proton would 
differ, since only for large y are G p p (y) ~ (1 — y) 5 jy and 
Gp(y) ~ (1 — y) 2 jy significantly different. In Ref. [15] we 
have shown that the introduction of the above motioned 
cuts drastically reduces the energy (y/s) dependence of 
the diffractive mass distributions leading, in particular, 
to the approximate \jM\ behaviour for all values of y/s 
from ISR to Tevatron energies. Here these cuts produce 
another type of scaling, which may be called "projectile 
scaling" or "projectile universality of the diffractive peak" 
and which means that for large enough xl the diffractive 
peak is the same for all projectiles. The corresponding \ d 



x For example, in the cosmic ray experiments it is usually 
assumed that K-kN = 1.5 K pp , which is traced to analysis of 
data like those in [69] performed in terms of the additive quark 
models (cf., [70]). 



functions will be the same for protons and pions in this 
region. The cross section appearing in the denominator 
of the moments will, in this case, be the same, i.e., a Fp . 

The only remaining difference between them, their dif- 
ferent gluonic distributions, is in this region cut off. This 
may be regarded as a prediction of the IGM. Experimen- 
tally this may be difficult to check since one would need 
a large number of points in large xl region of the leading 
particle spectrum. Data plotted in Fig. 18 neither prove 
nor disprove this conjecture. The discrepancy observed 
in the proton spectrum is only due to our choice of nor- 
malization of the diffractive and non-diffractive curves. 
The peak shapes are similar. 




p% + ^>pn + kad X . 
, i , , , i , , . i 
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FIG. 19. a) Comparison of our spectra F(xl) for 
leading pions with data from Ref. [62] in the reaction 
7r + + p — > 7r + + p + X. Solid and dashed lines correspond 
to the choices mo = 0.35 GeV and mo = 0.45 GeV respec- 
tively, b) the same as a) for the leading proton spectrum 
F(xf) measured in the same reaction. 
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The EHS/NA22 collaboration provided us with data 
on 7r + + p — > n + + p + X reactions. In particular they 
present the xf distributions of both leading particles, the 
pion and the proton. Their points for pions and protons 
are shown in Fig. 19a) and b) respectively. These points 
are presumably free from diffractive dissociation. The 
above mentioned asymmetry in pion and proton energy 
loss emerges clearly, the pions being much slower. The 
proton distribution peaks at xp — 0.6 — 0.8. Our curves 
(solid lines) reproduce with no free parameter this be- 
haviour and we obtain a good agreement with the pion 
spectrum. Proton data show an excess at large xp that 
we are not able to reproduce keeping the same values of 
parameters as before. 

The authors of Ref. [62] tried to fit their measured pro- 
ton spectrum with the FRITIOF code and could not ob- 
tain a good description of data. This indicates that these 
large xf points are a problem for standard multiparticle 
production models as well. In our case, if we change 
our parameter mo from the usual value mo = 0.35 GeV 
(solid line) to m = 0.45 GeV (dashed line) we can re- 
produce most of data points both for pions and protons 
as well. This is not a big change and indicates that the 
model would be able to accomodate this new experimen- 
tal information. Of course, a definite statement about 
the subject would require a global refitting procedure, 
which is not our main concern now. 



B. Leading particles in photon-proton collisions 

If, at high energies, the reactions p — p and it — p have 
the same characteristics and if VDM is a good hypothesis, 
then more about the energy flow in meson-p collisions can 
be learned at HERA. Indeed, as mentioned in [54], at 
the HERA electron-proton collider the bulk of the cross 
section corresponds to photoproduction, in which a beam 
electron is scattered through a very small angle and a 
quasi-real photon interacts with the proton. Using VDM, 
high energy photoproduction exhibits therefore similar 
characteristics to hadron-hadron interactions. 

Data taken by the ZEUS collaboration at HERA [63] 
show that the LP spectra measured in photoproducion 
and in DIS (where Q 2 > 4 GeV 2 ) are very similar, spe- 
cially in the large xl region. This suggests that, as 
pointed out in [72], the QCD hardness scale for parti- 
cle production in DIS gradually decreases from a (large) 
Q , which is relevant in the photon fragmentation region, 
to a soft scale in the proton fragmentation region, which 
is the one considered here. We can therefore expect a 
similarity of the inclusive spectra of the leading protons 
in high energy hadron-proton collisions, discussed above, 
and in virtual photon-proton collisions. In other words, 
we may say that the photon is neither resolving nor being 
resolved by the fast emerging protons. This implies that 
these reactions are dominated by some non-perturbative 
mechanism. This is confirmed by the failure of perturba- 



tive QCD [73], (implemented by the Monte Carlo codes 
ARIADNE and HERWIG) when applied to the proton 
fragmentation region. In Ref. [72] the LP spectra were 
studied in the context of meson and Pomeron exchanges. 
Here we use the vector meson dominance hypothesis and 
describe leading proton production in the same way as 
done for hadron-hadron collisions. The only change is 
that now we have p— p instead oip—p collisions. Whereas 
this may be generally true for photoproduction, it re- 
mains an approximation for DIS, valid in the large xl 
region. 

Assuming that VDM is correct, the incoming photon 
line can be replaced by solid line in Fig. 1. During the 
interaction the photon is converted into a hadronic state, 
called V, and then interacts with the incoming proton. 
At HERA only collisions V — p are relevant. The state 
V looses fraction x of its original momentum and gets 
excited carrying ajF = l- i fraction of the initial mo- 
mentum. The proton, which we call here the diffracted 
proton, looses only a fraction y of its momentum but oth- 
erwise remains intact. We assume here, for simplicity, 
that the vector meson is a p° and take G p (x) = G 7T (x) 
in eqs. (5) and (6). 

10 r ' 1 ■ ■ ■ 1 ' ' ' 1 ■ ■ ' 1 ' ■ ■ 1 

: . ZEUS 2002 

' 3 < Q 2 < 254 GeV 2 
45 < W < 225 GeV 

p 2 < 0.5 GeV 2 ' 

1 da 
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FIG. 20. Comparison between our calculation and the data 
on the leading proton spectrum measured at HERA by the 
ZEUS Collab. [23]. 

In Fig. 20 we present our spectrum of leading protons 
in 7p collisions. All parameters leading to the results 
in that figure arc the same as established before in our 
study of diffractive mass distributions in photon-proton 
collision at HERA. 



C. Leading J/ip production 

All produced particles come essentially from the glu- 
ons and quark-antiquark pairs already pre-existing in the 
projectile and target, or radiated during the collision. 
This qualitative picture takes different implementations 
in the many existing multiparticle production models. In 
the IGM, the produced particles (and consequently the 
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energy released in the secondaries and lost by the projec- 
tiles) come almost entirely from the pre-existing gluons 
in the incoming hadrons. This conjecture may be directly 
tested using a high energy, nearly gluonless hadronic pro- 
jectile. In this case, according to the IGM, inspite of 
the high energy involved, the production of secondaries 
would be suppressed in comparison to the production 
observed in reactions induced by ordinary hadrons. The 
energy would be mostly carried away by the projectile 
leading particle which would then be observed with a 
hard xp spectrum. This type of gluonless projectile is 
available in J/ip photoproduction, where the photon can 
be understood as a virtual cc pair which reacts with the 
proton and turns into the finally observed J /if). There 
are low energy data taken by the FTPS Collaboration 
[74] and high energy data from HERA [75]. 

We want to stress here the fact that the fair agreement 
with data observed in Figs. (18a) and (20) is possible 
only because the diffraction processes have been prop- 
erly incorporated in the calculations. In other words, the 
inclusion of a diffractive component turns out to be a 
decisive factor to get agreement with data. We can also 
describe reasonably well pionic and kaonic LP and the 
observed difference turns out to be due to their different 
gluonic distributions. 
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FIG. 21. Comparison of the IGM distribution F{z) with 
data of Ref. [75] with restricted acceptance v% > 1 (GeV/c) 2 
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three different values oip J/ *: 0.066 (dashed line), 0.033 (solid 
line) and 0.016 (dotted line). 



The crucial role played by the parameter p (see eq. (8)) 
representing the energy-momentum fraction of a given 
hadron allocated to gluons is best seen in Fig. 21 where 
we show the fit to data for leading J/ij) photoproduction 
[75] . The only parameter to which results are really sen- 
sitive is p = p J l^ which, as shown in Fig. 21, has to be 
very small, p J /^ = 0.033. This is what could be expected 
from the fact that charmonium is a non-relativistic sys- 
tem and almost all its mass comes from the quark masses 
leaving therefore only a small fraction, 



M 



j/i> 



2m c 



M 



0.033, 



(21) 



for gluons (here m c = 1.5 GeV and Mj^, — 3.1 GeV). 
Of course, the value of p J /^ required to give a very good 
fit of data might change either with another choice of 
m c or another choice of . However these changes 

might affect p J >^ by, at most, a factor two. This suggests 
that the momentum fraction carried by gluons in the J/tp 
is one order of magnitude smaller than that carried by 
gluons in light hadrons. 



IX. SUMMARY AND CONCLUSIONS 

We were able to fit an impressive amount of experimen- 
tal data, which had nothing in common except the fact 
that they always referred to the momentum (or rapidity) 
distribution of some observed particle or to the invariant 
mass distribution of a cluster of measured particles. We 
could fit these data starting from one single "generat- 
ing" function, xi x i y)i which depends almost only on the 
density and interaction cross section of the gluons inside 
hadrons. These are fundamental quantities in QCD and 
with our model we can test the existing results for them. 
More than just fitting, we did some predictions and one 
of them, the leading particle spectrum shown in Fig. 20, 
was confirmed by experiment. 

After all these works, we may ask ourselves what have 
we learned. We believe that we have constructed a sim- 
ple and consistent picture of energy flow in strong inter- 
actions, based on the assumption that energy loss and 
leading particle spectra are determined by many indepen- 
dent gluon-gluon collisions and valence quarks play a sec- 
ondary role. Consequently, energy flow will reflect the 
properties of the gluon distributions and cross sections 
in the colliding hadrons. This picture seems to be uni- 
versal, i.e., valid in many different contexts. However, in 
order to see this universality we have to be careful and 
use proper kinematical limits of the phase space for ev- 
ery reaction considered, as illustrated in Figs. 1 and 2. 
When this is done the sensitivity of energy flow to other 
(than gluon distributions and cross sections) aspects of 
the production process is only of secondary importance 
and needs special observables (which are sensitive to, for 
example, the quantum numbers of the detected particles) 
to be visible. But even then, the IGM is indispensable 
because it provides the important energy correlations be- 
tween different parts of the phase space. 

Our analysis shows also clearly that our model can 
be regarded as a useful reference point for all more so- 
phisticated approaches whereas, for hydrodynamical ap- 
proaches of multiparticlc production, it provides the ini- 
tial energy used for the further evolution and hadroniza- 
tion of the created systems. However, in order to comply 
with the recent developments of QCD concerning the low 
x gluonic content of hadrons [26,27] it must be accord- 
ingly updated. We plan to do this in the future. We also 
plan to account for the intrinsic fluctuations present in 
the hadronizing systems. In the usual statistical mod- 
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els this can be done by using the so called nonexten- 
sive statistics and, as was shown in [76], it can influence 
substantially some energy flow results, in particular the 
estimation of inelasticity K. 



APPENDIX A 



The main ideas 

The IGM can be summarized in the following way: 

(i) The two colliding hadrons are represented by va- 
lence quarks carrying their quantum numbers plus 
the accompanying clouds of gluons. 

(ii) In the course of a collision the gluonic clouds in- 
teract strongly depositing in the central region of 
the reaction fractions x and y of the initial energy- 
momenta of the respective projectiles in the form 
of a gluonic Central Fireball (CF). 

(Hi) The valence quarks get excited and form Leading 
Jets (LJ's) which decay and populate mainly the 
fragmentation regions of the reaction. 

The fraction of energy stored in the CF is therefore 
equal to K = ^fxy and its rapidity is Y = ^ In | . 

The CF consists of minifireballs (MF) formed from 
pairs of colliding gluons. In collisions at higher scales a 
MF is the same as a pair of minijets or jets. In the study 
of energy flow the details of fragmentation and hadron 
production arc not important. Most of the MF's will be 
in the central region and we assume that they coalesce 
forming the CF. The collisions leading to MF's occur at 
different energy scales given by Q\ = X{ yi s, where the 
index i labels a particular kinematic configuration where 
the gluon from the projectile has momentum Xi and the 
gluon from the target has yt . We have to choose the scale 
where we start to use perturbative QCD. Below this value 
we have to assume that we can still talk about individ- 
ual soft gluons and due to the short correlation length 
between them they still interact mostly pairwise. In this 
region we can no longer use the distribution functions ex- 
tracted from DIS nor the perturbative elementary cross 
sections. 

The central formula 



The central quantity in the IGM is the probability to 
form a CF carrying momentum fractions x and y of two 
colliding hadrons. It is defined as the sum over an unde- 
fined number n of MF's: 



x{x, y) = E ■ ■ ■ E 6 ^ ~ ni Xl 

x 5 [y - m 7/1 - • 
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n l y l 



(Al) 



The delta functions in the above formula garantee en- 
ergy momentum conservation and P(ni) is the probabil- 
ity to have collisions between gluons with Xi and j/j. 
The expression above is quite general. It becomes spe- 
cific when we define P(n,). The assumption of multiple 
parton-parton incoherent scattering (which is also used 
in Rcfs. [1,29-31,34]) implies a Poissonian distribution of 
the number of parton-parton collisions and thus P(rij) is 
given by: 



P{ni) = 



{n,i) ni exp(-ni) 



(A2) 



Inserting P(rij) in (Al) and using the following integral 
representations for the delta functions: 
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we can perform all summations and products arriving at: 

x(x, y) = -jtt-^ I dt du ex P W x + u v)\ 



x exp ■ 



{ni \e-^ tXi+uyi ^ - l] } | (A5) 



Taking now the continuum limit: 
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we obtain: 
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This function ui{x' ,y') is called the spectral function 
and represents the average number of gluon-gluon colli- 
sions as a function of x' e y'. It contains all the dynamical 
inputs of the model and has the form: 



ui(x',y') = 



Oggjx'y's) 

a(s) 



G{x')G{y') 



x 9 (x'y 1 - K 2 mtn ) , 



(A9) 



where G's denote the gluon distribution functions in the 
corresponding projectiles and a gg and a are the gluon- 
gluon and hadron-hadron cross sections, respectively. In 
the above expression x' and y' are the fractional momenta 
of two gluons coming from the projectile and from the 
target whereas K m i n = mo/ \fs, with mo being the mass 
of lightest produced state and y/s the total c.m.s. energy. 
m is a parameter of the model. 

The integral in the second line of eq. (A7) is dominated 
by the low x' and y' region. Considering the singular 
behavior of the G(x) distributions at the origin we make 
the following approximation: 

e -i(tx'W) _ i ~ + uy >) _ + uy 'f (A10) 

With this approximation it is possible to perform the 
integrations in (A7) and obtain the final expression for 
x(x,y) discussed in the main text: 
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where 



and 



D xy = (x 2 )(y 2 )-(xyf 



(x n y m )= f dxx n [ dyy m uj(x,y), (A12) 
Jo Jo 

Xo is a normalization factor defined by the condition: 

[ dx [ dy X (x,y)9(xy-K 2 lm ) = l (A13) 
Jo Jo 



The numerical inputs 

In order to evaluate the distribution (All) we need 
to choose the value of mo , the semihard scale pt min and 
define G(x) and a gg in both interaction regimes. We take 
PTmin = 2.3 GeV and m — 0.35 GeV. These are the 
two scales present in the model. The semihard gluon- 
gluon cross section is taken, at order a 2 , to be: 



Vqq{x, V, S) = K 



16p 2 T , 
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where 

H = 36T 
and 
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The parameter k is the one frequently used to incorporate 
higher corrections in a s and is 1.1 < k < 2.5 according 
to the choice of G(x), of the scale Q 2 and prmin- For 
PTmin = 2.3 GeF, k = 2.5. 

The coupling constant is given by: 



a s {Q 2 ) = 



12tt 



(33 - 27V/) In 



(A17) 



where A = 0.2 GeV and Nf — 3 is the number of active 
flavors. As usual in minijet physics we choose Q 2 — 
PTmin an d use the distributions G(x,Q 2 ) parametrized 
in literature. 

When the invariant energy of the gluon pair s is the 
interval < s = xys < ^p\ min we are outside the 
perturbative domain. Parton-parton cross sections in the 
non-perturbative regime have been parametrized in [77] 
leading to a successful quark-gluon model for elastic and 
diffractive scattering. Recently these non-perturbative 
cross sections have been calculated in the stochastic vac- 
uum model [78]. The obtained cross sections are func- 
tions of the gluon condensate and of the gluon field corre- 
lation length, both quantities extracted from lattice QCD 
calculations. In order to keep our treatment simple we 
adopt the older parametrization for the gluon-gluon cross 
section used in [77]: 



K q {x, y, s) 



xys 



(A18) 



where a is a parameter of the model [11]. 

The main distributions 

Given \(x, y) we can immediately write the inelasticity 
distribution, its complementary distribution, the leading 
jet momentum spectrum and the CF rapidity distribu- 
tion: 



X(K) 



f dx f 
Jo Jo 



dyx(x,y)S(V x ~y- K ) 



x 9 (xy - K 2 mm ) 



(A19) 



F(x L ) = (l-a)[ dx X (nd) (x-y = l-x L ) + 

+ a i I dx X {d} {x;y = l-x L ) (A20) 

j=l,2 Jx ™- 
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^ Y ) = J o dx dyx(x,y)s(^ 
x 0(xy-K^i n ) 



In — 



Y 



(A21) 



where a = a\+a2 is the total fraction of single diffractive 
(d) events (from the upper and lower legs in Fig. 1, 
respectively) and where 



Max 



(M LP + m )f 



(1 - xl)s 



(A22) 



The mass spectra for Single Diffractive processes are 
given by: 



1 da 



adM x 



dN 
dM\ 



dx I dy X (x,y) 
o Jo 



6 (M x - S y) (xy - K, 



2 ■ ) 

rntn } 



(A23) 



Substituting now eq. (All) into eq. (A23) we arrive 
at the following simple expression for the diffractive mass 
distribution: 



' ;V = ±F(M x ,s)H(M x ,s) 



dM x 



where 

F(M 2 x ,s) =cxp 

and 

H(M x ,s)- X] 



(x 2 ) 
2D 



xy 



(y) 



(A24) 



(A25) 



f,*"»{-2k; z }- ,A26) 



Z = \tf)(x - {x) f - 2(xy)(x - {x))(M%/s - {y))] 

We first keep only the most singular parts of the glu- 
onic distributions used (i.e., G(x) ~ 1/x) and collect all 
other factors in eq. (A9) in a single parameter c. Assum- 
ing that the ratio of the cross sections a< f^ does not 
depend on x and y and neglecting all terms of the order 

2 2 

of ^ and -Sf- , we arrive at the following expressions for 

m x 

the moments calculated in eq. (A12): 



(x) = 2 (a; 2 ) ~ c In 

s , 2 , Ml , Ml 



(A27) 
(A28) 



, , {Ml ml , Ml \ , A . 

(xy) ~ c ^ In — ^ A29 

V s s m o / 

Notice that in all cases of interest (x y) is much smaller 
than other moments (by a factor hi^r, at least). It 
means that D xy ~ (x 2 }(y 2 ) and consequently 



F(A4,s) ~ exp 



2(y 2 ) 



2n 



1 — c In 



Jit 



and 

ff(Mf, S ) 



exp 



Xo 



c In 



(A30) 



2^^ 



dx exp 



(x-(x» 2 



2(x 2 ) 



v/^ 5 ) 1 

const - !L -= = const ■ 



~ const • 



M x Jc\n^kr 



(A31) 



leading to 



<i/V 1 , 9 N , 9 const 
r * -tf(M 2 , S )F(M 2 , S )~ x 



x 



exp 



In 



Mi 

v m o / 

c in — <f- 



(A32) 



The mass spectra for Double Pomeron Exchange pro- 
cesses are given by: 



1 da 



dN 



dx dy X (x,y) 
Jo Jo 



a dM X Y dM X Y 

x 5 (M XY - ^xyl) 9 {xy - K 2 mm ) (A33) 



APPENDIX B 



The "kinematical" Pomeron 

The Pomeron for us is just a collection of gluons which 
belong to the diffracted proton or antiproton. These 
gluons behave like all other ordinary gluons in the pro- 
ton and have therefore the same momentum distribution. 
The only difference is the momentum sum rule, which for 
the gluons in IP is: 



/' 

Jo 



dx' x' Gp(x') = p d 



(Bl) 



where = 0.05 (see [15]) instead of p — 0.5, which holds 
for the entire gluon population in the proton. 

In order to make contact with the analysis performed 
by HERA experimental groups we consider two possible 
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momentum distributions for the gluons inside P. A hard 
one: 

fg/M=a h {l-0) (B2) 

and a "super-hard" (or "leading gluon" ) one: 

f$ P (P) = a ah f?{\ - /3) - 3 (B3) 

where j3 is the momentum fraction of the Pomeron carried 
by the gluons and the superscripts h and sh denote hard 
and superhard respectively. The constants ah and a s h 
will be fixed by the sum rule (Bl). In the past, following 
the same formalism, we have also considered a soft gluon 
distribution for the Pomeron of the type 

(1 - 0) 5 

f S a /M= Q -^ EL ( B4 ) 

but we found that this "soft Pomeron" distribution was 
incompatible with the single diffractive mass spectra 
measured at HERA [79]. This Pomeron profile was also 
ruled out by other types of observables, as concluded in 
Rcfs. [80] and [54]. 

We use the Donnachie-Landshoff Pomeron flux factor, 
which, after the integration in the t variable, is approxi- 
mately given by [58]: 

f P/p (x P ) ~ Cx^ 2 "' ~ C — (B5) 

Xp 

where xp is the fraction of the proton momentum carried 
by the Pomeron and the normalization constant C fixed 
also with the help of (Bl). Noticing that (3 = the 
distribution Gp(y) needed in eqs. (5) and (6) is then 
given by the convolution: 

G s p h ' sh (y) = f^r- f P / P M f S g/ p h (—) (B6) 

Jy X]P X]P 

We use also the "diffractive gluon distribution" given 
by: 

(1 - v) 5 

G F (y)=a [ - y ±- (B7) 

V 

where a is fixed by the sum rule. 
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